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Abstract 

We present a detailed study of the decays of the neutral Higgs bosons 
and A^ within the Minimal Supersymmetric Standard Model. We find that the 
supersymmetric modes it and bb can dominate the and A^ decays in a wide 
range of the model parameters due to large Yukawa couplings and mixings of t 
and b . Compared to the conventional modes tt and bb, these modes have very 
distinctive signatures. This could have a decisive impact on the neutral Higgs 
boson searches at future colliders. 



The existence of two or more Higgs bosons would be a clear indication that 
the Standard Model must be extended. For example, the Minimal Supersymmetric 
Standard Model (MSSM) with two Higgs doublets predicts the existence of five 
physical Higgs bosons hP , H^, A^, and H"^ 0? §]• In searching for the Higgs bosons 
it is necessary to study not only the production mechanism, but also all their pos- 
sible decay modes. Their decays to supersymmetric (SUSY) particles could be very 
important if they are kinematically allowed. This is indeed the case for the charged 
Higgs boson H~^. If all SUSY particles are very heavy, the decays dominantly 
into tb ; the decays —>■ t^u and/or W^h^ are dominant below the tb 

threshold p|, §]. In refs.|^ H all decay modes of including the SUSY-particle 
modes were studied in the the case that the SUSY-particles are relatively light ; it 
was shown that the SUSY decay mode tb can be dominant in a large region of 

the MSSM parameter space due to large t and b quark Yukawa couplings and large t- 
and 6-mixings, and that this could have a decisive impact on searches at future 
colliders. In this paper we extend this study to the decay of the neutral Higgs bosons 
i/o and 

The lighter stop ti can be much lighter than the other squarks and even lighter 
than the t quark due to large ti-tR mixing being proportional to the large top Yukawa 
coupling ht and the t-mixing parameters At and fi Similarly, the lighter sbottom 
bi can also be much lighter than the other squarks In the case of large t- and 
b- mixings one also expects the couplings of H'^tt, A^tt, H^bb and A^bb to be large. 
Here we show explicitly that the modes tt and bb can indeed dominate the and 
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A" decays in a wide range of the MSSM parameters. 

First we summarize current experimental limits on the squark masses. The D0 
group at FNAL [H, |[ obtained mass bounds for squarks and gluinos from the negative 
search for them. However, they obtained no mass bound for the mass-degenerate five 
flavors of squarks (g ) (other than the stop (t )) within the MSSM in case the lightest 
neutralino Xi (assumed to be the lightest SUSY particle (LSP)) is heavier than ~ 
70GeV ; this is mainly due to the massive LSP effect. On the other hand, the LEP-I 
experiments set the general bound on the q mass mq ^45GeV |TI| . Experiments 
dedicated to t search at LEP1.5 ||TT| and Tevatron |]T^ have set some limits on the 
lighter stop (ti) mass for relatively light Xi? e-g-; typically m^^^SSGeV (for no t - 
mixing) and m^^^lOOGeV ||l2l, respectively. However, they set no limit on m^^ 



for m^o^SOGeV 0,|T2|. 

In the MSSM the properties of the charginos xf (^ = 1,2) and neutralinos 
(j = 1, ■ ■ ■ , 4) are completely determined by the parameters M, and tan j3 = V2/V1, 
assuming M' = (5/3) tan^ 9wM. Here M (M') is the SU(2) (U(l)) gaugino mass, /i is 
the higgsino mass parameter, and f 1 (f 2) is the vacuum expectation value of the Higgs 
Hi {H2) HT]. Here < m-+ and m^o < ■ ■ ■ < m^o. To specify the squark sector 
additional (soft SUSY breaking) parameters Mq, Mq, M^, (for each generation) and 
A (for each flavor) are necessary. The mass matrix for stops in the base of [t^, tn) 
reads 0, |] 

Ml=(< (1) 
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with 



mf^ = M| + + m| cos 2p{lf - tt sin^ dw), (2) 
m^^ = M|- + + m| cos 2/?ef sin^ ^vy, and (3) 
a^mt = mt(ylf - /icot/3). (4) 

For the 6 system analogous formulae hold but with M?. replaced by in eq.(3), 
and instead of eq.(4), abmf, = mfe(Af, — yUtan/5). mixing may be important for 

large Ay,, fi and tan/5. The slepton sector is fixed by adding (soft SUSY breaking) 
parameters and M^. 

The masses and couplings of the Higgs bosons H^, H^, and ^4°, including 
leading Yukawa corrections, are fixed by m^o, tan/3, rrit, Mq, Mfy, M^, At, Af, and /i. 
if'' and A'' are the heavier (lighter) CP-even and CP-odd neutral Higgs bosons, 
respectively. For the Yukawa corrections to the and masses and their mixing 



angle a we use the formulae of ref.|13]. We include also the potentially large Yukawa 



corrections to the triple Higgs vertices of H^h^h^ , H^A^A^ and h^A^A^ using the 



formulae of ref.in 



In the following, we take for simplicity Mq = Mq = Mf) (for the third gener- 
ation). Ml = M^ = Mq, and At = Af, = A^ = A. Thus we have as free parameters 
rriAO, M, fj,, tan/3, Mq and A. 

We calculate the widths of all possibly important modes of and A^ decays: 
(i) H'^ tt, hi, cc, T-T+, W+W-, Z^Z^, h^h^, A^A^, W^H^, Z^A^, Uj, bibj, £-£+, 
UiUi {i = e,ii,T), xtXj, XkXh and (ii) A^ ti, hi, cc, t't^, Z^h°, 1^2, hh, 6162, 



t^iT^2j ^2^1 ^ xtxjj Xk'Xj- Formulae for these widths are found in ref. 0. In 
principle, also the decays QaQa il = ""^ c, s and a = L,R) could contribute 

via their gauge couplings. As the squarks of the first two generations are supposed 
to be heavy, these decays will be strongly phase-space suppressed. Even if they were 
kinematically allowed, they would have a rate at most comparable to that of if ° — 

and ui^i (see fig. 2 below). We neglect loop induced decay modes (such as 



H — > gg and 77) and three-body decay modes |]15 . 

In order not to vary many parameters, in the following we fix rrit = ISOGeV 
and /i = SOOGeV, and take the values of M and tan/? such that m^o ~ 70GeV 



for which the D0 bounds on and m^-^ p, 12] and the LEP1.5 bound on m^-^ |11 
disappear. 

In fig.l the contour lines for the branching ratios B{H —>■ tt, bb) = 
Ei,j=i,2iBiH Uij)+B{H bibj)) {H = H° or A°) are plotted in the A-Mq plane 
for (mAo(GeV), M(GeV), tan/5) = (450, 160, 2) (a), (500, 160, 2)(b), (450, 146, 12) (c) 
and (500, 146, 12)(d), for which m-+(GeV) = 128, 128, 131 and 131, respectively. In 
the plots we have required iii^ ^xl (~ 70GeV). In the allowed regions of fig.l 
we find mijo ~ 771^0 and m/jo > 55GeV, which satisfies the LEP limit m^o > 44GeV 



| T6[] . We see that the branching ratios B{H^ tt,bb) and 5(^4° tt,bb) can be 



larger than 50% in a sizable region. In this region, these SUSY decay modes dominate 
over the conventional modes. 

In fig. 2 we show the m^o dependence (in the rriAO range for rriho > 50GeV) 
of the important branching ratios of and ^4° decays for (MQ(GeV), y4(GeV), 
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M(GeV), tan/3) = (102, 325, 160, 2) (a, b), (218, 405, 145, 30) (c), and (145, 274, 
146, 12) (d). In these three cases we have {rrif^, rrif^, m^^, mi^, ^^+) = (100, 270, 
100, 115, 128) GeV (a, b), (80, 386, 80, 304, 132) GeV (c), and (80, 310, 80, 199, 131) 
GeV (d) . We see that in these cases the sum of the tt and bb modes dominates the 
and ^4*^ decays in a wide range of m^o. Here note that ttiho — m^o in the m^o 
range shown here, and that the A° does not couple to fjij, bibi and f^ff {i — 1,2). 

As for h° decay, we have found that the decay hP — > titi is kinematically 
allowed only in a very limited region of the MSSM parameter space. 

The reason for the dominance of the tt and bb modes in the decay is as 
follows : The modes tt, bb, tt and bb (whose couphngs to are essentially ~ ht sin a, 
~ hbcosa, ^ {A — 11 cot a)ht sina and ~ (^4 — iJ,teiRa)hbCosa, respectively) can be 
strongly enhanced relative to the other modes due to the large Yukawa couplings hf^. 
In addition, the tt and bb modes can dominate over the tt and bb modes, respectively, 
in the case the ^-mixing parameters A and /i are large. Furthermore, in this case ii 
and bi tend to be light due to a large mass-splitting. Here note that the effects of the 
bottom Yukawa coupling ht, and b- mixing play a very important role for large tan /3 
(see fig. 2c). The reason for the dominance of the ii and bb modes in the A° decay is 
similar to that in the decay. 

Quite generally, B{H — > ii,bb) depends on the parameters m^o, Mq, A, /i, 
tan P and more weakly on M. For a given m^o the strongest dependence is that on 
Mq to which ruf and are sensitive (see fig.l). B{H — > ii,bb) can be quite large 
in a substantial part of the parameter region kinematically allowed for the ii and bb 
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modes. We find that the dominance of the tt and bb modes is fairly insensitive to the 
assumption Mq = Mfj = Mfj = Mi = M^. As seen in fig.l the dependence on A is 
also strong. Concerning the assumption At = Ai, = A^, we have found no significant 
change of B{H —>■ tt,bb) {H = H^,A^) as compared to fig. 2, when we take Ab^r/At 
= ±0.5, ±1, ±2 keeping A^ = A. As \n\ increases, the dominance of the it and bb 
modes becomes more pronounced because of the increase of the {H^, A^) couphngs 
to tt and bb. We also find that B{H tt, bb)) {H = H^, is nearly invariant under 
(/i. A) ^ i-fx-A). 

Here we have not included QCD radiative corrections to the hadronic modes. 



It is shown in ref.[17| that both the standard QCD correction (due to gluon-quark 
loop) and SUSY QCD correction (due to gluino-squark loop) to the widths of H^, A^ 
—>■ tt, bb can be large, but that the two corrections can partly or even totally cancel 
each other. The QCD corrections to H^, A^ — > tt, bb are not known. It is shown in 
ref. that the QCD corrections to the closely related process —>■ tb are significant, 
but that they do not invalidate the tree- level conclusion of ref . |^ on the dominance of 
the tb mode in a wide parameter region. As the structures of the couplings of and 
A'^ to it and bb (ti and hb) are similar to those of to tb {tb), we can expect that 
the QCD corrections would not invalidate our tree-level conclusion on the dominance 
of the tt and bb modes in a large parameter region. 

As for the signatures of the and A^ decays, typical tt and bb signals are 
shown in table 1. They have to be compared with the conventional ti and bb signals, 
respectively, (if°, A^) ^ ti ^ {W+b){W-h) {ff')b{ff')b (i.e., 6 jets (j's), 4j's + 1 
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isolated charged lepton (i^) + missing energy-momentum (^), or 2j's + + />) and 

y4°) bb (i.e., 2j's). Note that B{ii cXi) ^ 1 if m^-^ < m-+ ^- - and m^o < 
< m^o + rrit^w (in cases (a), (b) and (d) of table 1), and B{ii cx?) — otherwise (in 
case (c)) [|18|. (In principle there is a region mb + mw + m^o < mi^ < mb + m^+ where 
also the decay ii — * bW^Xi plays a role However, this decay does not occur in 
the cases considered here.) As seen in table 1, the tt (bb) signals have general features 
which distinguish them from the tt (bb) signals : (i) more p due to the emission of 
two LSP's and hence less energy-momentum of jets and the isolated-charged-lepton 
in case of a short decay chain, or (ii) more jets and/or more isolated-charged-leptons 
in case of a larger decay chain. Moreover, depending on the values of the MSSM 
parameters, the tt and bb signals could have the following remarkable features : (i) 
production of same-sign dileptons (e.g., in case (g)), which could yield same- 

sign isolated dilepton events such as e~^e~ — * H^Z^ (£+/+ or £^£'^) + j's + ^ ; (ii) 
less bottom-jet activity (e.g., in cases (a), (b), (d) with (/i°, Z°) (£+£^ or uu) ) or 
more bottom-jet activity (e.g., in cases (c), (d), (f) with (/i°, Z^*^) bb) ; and (iii) 
emission of a real or (e.g., in cases (b), (c), (d), (f)). The identification of the 
sign of charged leptons and the tagging of b- and c-quark jets, h'^, Z^ and would 
be very useful in discriminating the tt and bb signals from the tt and 66 signals as 
well as in suppressing the background. 

The suitable places for and search would be e^e~ colliders ||2^ and 
hadron supercolliders 0]. If the tt and bb modes dominate the if" and decays, 
it decisively influences the signatures of and . For example, it can strongly 
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suppress the conventional discovery modes of and at the hadron supercolhder 
0], such as the "gold-plated" —>■ Z^Z^ mode (viable for small tan/5 and 2mz 
<, rriHQ <, 2mt) (see fig. 2a) and t^t^ mode (viable for large tan/?) (see 

figs. 2c, d). Clearly it would be necessary to perform a detailed Monte-Carlo study 
to separate the signals from the background. Such a study is, however, beyond the 
scope of this article. 

We have shown that the SUSY modes tt and bb can dominate the and A^ 
decays in a large allowed region of the MSSM parameter space due to large t and b 
quark Yukawa couplings and large t - and b -mixings. The tt and bb modes have very 
distinctive signatures as compared to the conventional modes tt, bb and r~r+. This 
could decisively influence the and A^ search at future e^e^ and hadron colliders. 

While preparing this manuscript a paper by Djouadi et al. |21| appeared deal- 
ing with a similar subject. They also point out the possible importance of the squark 
pair modes in and A^ decays. They, however, studied them in a strongly con- 
strained supergravity model in contrast to our work, which is in the general framework 
of the MSSM. 
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Figure Captions 



Fig.l Contour lines of B{H^ it, bb) (a, c) and B{A'^ it, bb) (b, d) in the A - Mq 
plane for (mAo(GeV), M(GeV), //(GeV), tan/3) = (450, 160, 300, 2)(a), (500, 
160, 300, 2)(b), (450, 146, 300, 12) (c) and (500, 146, 300, 12) (d). The contour 
oi B — (dashed hne) represents the kinematical boundary ; the hues (iii) , (iv) 
and (v) are the contour lines of 2m£^ = itiho (~ 450 GeV), 2m^^ = m^/o (~ 450 
GeV) and mg^ + m^^ = m^o (= 500 GeV), respectively. The shaded area is 
excluded by the requirement m^^ iii ^ ^xl (— GeV) ; the lines (i) and (ii) 
are the contour lines of m^^ and mg^ = 70 GeV, respectively, above which m^^ 
and > 70 GeV. 

Fig. 2 The m^o dependence of important branching ratios of the and decays 
for (MQ(GeV), ^(GeV), M(GeV), /x(GeV), tan/?) = (102, 325, 160, 300, 2) 
(a, b), (218, 405, 145, 300, 30) (c) and (145, 274, 146, 300, 12) (d). The sum 
over all mass eigenstates and/or flavors is taken for B{ii), B(bb), B(x'^x~): 
S(x°x°), ^(f+f") and B{i>u). The shoulders in the curves for the ii and bb 
modes correspond to opening of the ^2^2 and 62^2 channels, respectively. 
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Table 1: The typical tt and bb signals of the A°) decays in comparison to the 

conventional tt and bb signals, respectively. The decay chains (a), (b) and (d) are 
possible if m^^ < m-+ ^ - and m^o < m^^ < m^o + nit^w- i, -^^*'' and / denote 
missing energy-momentum, jet, isolated charged lepton, real (or virtual) Z° boson, 
and [q, £~ , v), respectively. A° does not couple to qiq^ and q'2<?2- 



Typical decay chains 



Signatures 



(a) H° 



ii 



2j's + i> 

(more f>, less h activity) 



2j's + i+i- + i, 
(Z° or h° emission ; 
less h activity if Z° or h° 
£+£- or lyjy) 



ii ^ ex? 



(b) R°,A° 



t'j 



t,Z° ^ (cxWf) 
or 

iih° (cx°i){hh or T+T-) 



4:j's + £+£- + ^, 
2j's + £+£-£+£- + i, 
2j's + £+£- + ^, 

or h° emission ; 
more h activity if Z° or h° 
hh ; trilepton) 



(c) E°,A° 



{ h ^ hxt ^ Hff'xl) 

i,z° ^{bxTXff) 
^ Hff'xDiff) 

or 

iih° (hxi)(hh or r+r") 
-^Hff'X°i)(hh or T+T-) 



6j's + ^, ij's + £+£- +i>, 
Aj's + £+ Aj's + ^, 
2j's + £+£-£+ + i, 
2j's + £+ 

(Z° or h° emission ; 
more h activity if Z° or h° 
bb ; less b activity if Z° 
or h° £~^£~ or vi? ; trilepton) 



(d) H"^ 



i2 ^ tx° ^ {bW+)x° ^ b{ff')x° 



iih° icxi)ibb or r+r") 



(e) R°,A° 




hxl 
hx'i 



2j's + ^ 
(more jd) 



(f) R°,A° 



hi 



bx 



ixt 



b{z(*)x\) ^ b{ff)x\ 

or 

b{h°xX)^b{bb or T+T-)x\ 

■{hw-){fTx\) 
Hfhff'xl 



8j's + ^, 6fs + £+£- 
Qi's + £^ +i>, 6j's+^, 
Ai's + £+£-£^ 

Aj's + £+£- + i), Aj's + £^ +i>, 
2j's + £+£-£+£- + ^, 
2j's + £+£- + i, 
(Z° or h° emission ; 
more b activity if h'^ or Z^*) 
bb ; trilepton) 



6i's + £±£± 

Ai's + £+£-£^ 

2j' s + £+£-£+ £- + j) 

(same-sign dilepton ; trilepton) 



(g) HO 



h ^ ^xr ^ (w+)(//'x?) ^ Hff')ff'xl 



&2 ^ ^ (W-)(//'X?) ^ b{ff')ff'x° 
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